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ABSTRACT 

We report on observations of giant pulses from the Crab pulsar performed simultaneously 
with the Parkes radio telescope and the incoherent combination of the Murchison Widefield 
Array (MWA) antenna tiles. The observations were performed over a duration of approximately 
one hour at a center frequency of 1382 MHz with 340 MHz bandwidth at Parkes, and at a center 
frequency of 193 MHz with 15 MHz bandwidth at the MWA. Our analysis has led to the detection 
of 55 giant pulses at the MWA and 2075 at Parkes above a threshold of 3.5cr and 6.5cr respectively. 
We detected 51% of the MWA giant pulses at the Parkes radio telescope, with spectral indices 
in the range of—3.6 > a > —4.9 (S^ oc u a ). We present a Monte Carlo analysis supporting 
the conjecture that the giant pulse emission in the Crab is intrinsically broadband, the less than 
100% correlation being due to the relative sensitivities of the two instruments and the width of 
the spectral index distribution. Our observations are consistent with the hypothesis that the 
spectral index of giant pulses is drawn from normal distribution of standard deviation 0.6, but 
with a mean that displays an evolution with frequency from -3.00 at 1382 MHz, to -2.85 at 192 
MHz. 

Subject headings: ISM: individual (Crab Nebula) - ISM: structure - pulsars: general - 
pulsars: individual (Crab pulsar) - scattering 
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1. INTRODUCTION 


The Crab pulsar (PSR B0531+21) w as first re¬ 


veale d by its exceptionally bright pulses ( Staelin fc Reifenstojgio 


1968ll . These so-called “giant pulses” are short du¬ 
ration (ranging from a few ns to a few ps) radio 
bursts, occurring only at the main-pulse (MP) 
or th e inter-pulse (IP) phases of the pulsar rota¬ 
tion (lHa nkins et^h 2003^ Hankinsfc bilekl [2007t 


Bhat et all 2008t Karuppusamv et all 2012h . The 

giant pulse amplitudes and energies typically ex¬ 
ceed those of regular pulses by two to four or¬ 
ders of magnitude. While the pulse energy dis¬ 
tribution of normal pulses fol lows an exponenti al 
or log-normal distribution (iBurke-Spolaor et al. 
20li), giant pulses instead we re originally ob- 
served to obey a power-law (LW£vle_feGower 


1972; iHesse fc Wielebinskil Il974t lR.itchingslll976t ) 
suggesting that they are created by a different 
mechanism. Howe ver, subsequent observations 
( Cordes et al. [2004 1 indicate that this power-law 
interpretation is an oversimplification only valid 
over small separations in radio frequency. 

The intrinsic properties of giant pulses are often 
obfuscated by interstellar scintillation and scatter¬ 
ing. For instance, observations at low frequencies 
are significantly affected by pulse broadening that 
arises from multi -path scattering in the interstel¬ 
lar m e dium (ISM; Bhat et al. 2007tlStappers et aT 


20111 ; lEllingson et al. I |2013j). Other than the 
well-known Crab, giant pulse emission has also 
been observed from the LMC pulsar PSR B0540- 
69 (J0540-6919), and a number of millisecond 
pulsars including PSR B1937+21 (01939+2134) 


(IWolszczan et al 


g rsti 

1. 1981 


__ Romani & Johnston " 2001 

Knight et alJ 120051) . Large amplitude pulses 
of millisecond widths have a lso been observed 
in so me long-period pulsars ( Ershov fo Kuzminl 
20051 . 

In the case of the Crab, giant-pulse emission 
is seen across the full range of the electromag¬ 
netic spectrum, from radio to gamma-rays. How¬ 
ever, the physics governing the e mission mech- 


anism is still not well understood (Cordes et al 


2004 ; Bilous et al.ll20li ). Inferred brightness tem¬ 
peratures are of the order of 10 30-32 K even when 
observations are affected by scattering or instru- 
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mental broadening, reaching up to ~10 41 K in 
observations made at ultra-high t i me resolutions 
(be, sub-ns) (lHankins et al. 20031: Kostvuk et al 


Cordes et al. 


Bhat et all 20081). 


20041 lHankins fc Eilekl l2007t 


The short duration of each giant pulse implies 
broadband e mission as infe rred from Heisenberg- 
Gabor limit ( Gabor 19461) . Simultaneous obser¬ 


vations of the Crab pulsar at widely separated 
frequencies have been empl oyed to investigate the 
validity of this assumption. Cornelia et al . ( 19691) 


first reported 50% correlation of giant pulses de¬ 
tected at 74 and 111 MHz. Different degrees of 
correlation have since been observed from as low 


as 3% to as high as 70% 

Sallmen et al. 19991 

Podov et al. ; 20061: Bhat et al 

20081: Mickaligeret al. 

2012). 


In order to investigate how correlated giant 
pulses are across wide frequency separations, we 
have performed an experiment using the Murchi¬ 
son Widefielcl Array (MWA) and the Parkes radio 
telescope. We report observations performed si¬ 
multaneously with the MWA, at a frequency of 
193 MHz (15 MHz bandwidth) and the Parkes ra¬ 
dio telescope, at a frequency of 1382 MHz (340 
MHz bandwidth). The MWA is a low frequency 
(80-300 MHz) interferometric array of 128 aper¬ 
ture array tiles designed for the detection of the 
redshifted neutral hydrogen 21-cm signal from the 
Epoch of Reionization (EoR), Galactic and ex- 
tragalactic surveys, Solar, heliospheric and iono¬ 
spheric studies, as well as transient and pulsar 
studies (Lonsdale et al. 112001): [Bowma n et al.j2013t 


Tingav et al.l [20130 . The MWA entered into its 
full operational phase in mid 2013. 

In this article, we independently determine 
the time-of-arrival of pulses for each instrument. 
Twenty three coincident giant pulses were de¬ 
tected. We compare our results to those in the 
literature, determine an intrinsic giant pulse spec¬ 
tral index distribution via a Monte Carlo analysis, 
and discuss the constraints we can place on giant 
pulse emission mechanisms. 

2. OBSERVATIONS 
2.1. Parkes observations 

The observations presented here were made on 
23 September 2013 with the 20-cm multibeam re¬ 
ceiver system on the Parkes radio telescope. These 
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data were recorded with a digital backend in use 
for the ongoing High Time Resolution Universe 
survey, the Berkeley - Park es Swinburne Recorder 
1BPSR: [Keith et al.ll2Q10h . from UTC 20 h 12 m 50 s 
to UTC 21 h 02 m 10 s , with a break of 245 s to al¬ 
low for system calibration (from UTC 20 h 39 m 05 s 
to UTC 20 h 43 m 10 s ). The data were recorded 
over a bandwidth of 340 MHz centered at 1382 
MHz at each of the two linear polarisations. The 
polyphase filters implemented on the FPGA in 
the BPSR channelize the digitized data into 1024 
channels with a spectral resolution of 390 kHz. 
The individual channels were subsequently de¬ 
tected and summed over both polarizations and 
25 samples to yield a 64 ps resolution time series. 
See Table [T] for a summary of the observation pa¬ 
rameters. 


2.2. MWA observations 

The MWA observations presented here were 
recorded from UTC 20 h ll m 08 s to UTC 20 h 55 m 06 s , 
using the MWA Voltage Capture System (VCS). 
The details of the VCS are presented in Tremblay 
et al. (2015). 


Table 1: Observation parameters. 


Parameters 

MWA 

Parkes 

T sys (K) 

268 

23 

Gain (KJy 1 ) 

0.012 a 

0.74 

HPBW 

24° 

14' 

Center Frequency (MHz) 

192.64 

1382 

Bandwidth (MHz) 

15.36 

340 

b Time Resolution (ps) 

400 

128 

Number of detected pulses 

55 

2075 

Average pulse rate 



(pulses/minute) 

1.25 

42.35 


a The net Gain is ~0.042 K Jy _1 for zenith 
pointing. 

b Resolution of the timeseries giant pulse search 
was performed on. 


We used the VCS to record complex voltages 
in 10 kHz channels at the Nyquist-Shannon rate, 
resulting in one complex voltage sample every 100 
ps. The 10 kHz channels were arranged in 12 
groups of 128 channels (each group representing 
a 1.28 MHz sub-band in the MWA system). From 
each 1.28 MHz sub-band, the central 88 channels 


were recorded owing to a limitation in our data 
recording during the commissioning of this observ¬ 
ing mode. Voltages for both linear polarizations 
were recorded for 15.36 MHz (half the maximum 
bandwidth of the MWA system) centered at 192.64 
MHz for all 128 tiles. See Table [l] for a summary 
of the observation parameters. 


3. DATA ANALYSIS AND RESULTS 


3.1. Parkes data analysis 


The data were searched in real-time for gi¬ 
ant pulses using the Heimdall single pulse pro¬ 
cessing software!^, following procedures s imilar to 
those described in Burke-Suolaor et all (2011). 
The giant pulse topocentric arrival times were 
recorded, along with other properties including 
signal to noise ratio, dispersion measure (DM), 
and matched filter width (128 ps). A total of 2075 
giant pulses above a signal to noise ratio (S/N) of 
6.5 were detected. This detection threshold is a 
default setting by the instrument’s hardware due 
to the RFI environment. 


Assuming the noise is Gaussian, we can esti¬ 
mate the number of false positives above the de¬ 
tection threshold using the Gaussian distribution 
function. The probability, P n that we obtain an 
event above a threshold of na is 


P n (x > na) 




P(x)dx = -erfc 


n 

.71 


(i) 


where a is the rms fluctuations in the noise, /i 
is the mean noise level, and erfc is the comple¬ 
mentary error function. The number of false pos¬ 
itives A, above a threshold of na is the product 
of the probability P ni of obtaining a false posi¬ 
tive above this threshold from one sample, and 
the total number of samples. Above 6.5cr, we es¬ 
timated the number of false positives to be signif¬ 
icantly less than one. Throughout the Parkes ob¬ 
servations, we detected 3 events above 6.5er in the 
off-pulse region. Which is higher than predicted, 
and probably indicative of radio frequency inter¬ 
ference. However the on-pulse region, to which we 
restrict our detections, is less than I percent of 
this extent. We therefore do not expect any false 
positive detections in our data set. 


1 http: //sourceforge.net/projects/heimdall-astro 
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The flux density, S of each giant pulse was cal¬ 
culated using the radiometer equation 


S = 


(S/N)S s 

Tlpti/ 


tAv 


( 2 ) 


where S sys is the system equivalent flux density, 
n p = 2 is the number of polarisations, ti nt = 128 
ps is integration time, and Av = 340 MHz is the 
observing bandwidth. Given that the size of the 
radio telescope beam (See Table [TJ is only about 
three times larger than the characteristic size of 
the Crab Nebula (~5.5'), the S sys is dominated by 
the flux contribu t ion fr om the Nebula. Following 
Bietcnholz ct al. ( 1997 b the flux contribution is 
Son = 955i' -0,2 ' Jy. The total system equivalent 
flux density is therefore 


Ssys — S syso A ScNi 


(3) 


where S syso = T syso /G is the system equivalent 
flux density in the absence of the Crab Nebula, 
T syso is the system temperature in the absence of 
the Nebula, and G is the gain of the telescope. 
This translates to a S sys of 906 Jy. 

Although pulsars are in general stable rota¬ 
tors, the Crab pulsar is not. It is a young pul¬ 
sar and subject to abrupt changes in rotation rate 
( Lvnc ct al.lll993r) . Therefore a precise ephemeris, 
valid for the day of observation, had to be obtained 
in order to measure the rotation rate of the pul¬ 
sar. In this case we used the time of arrival of the 
giant pulses at Parkes, all of which are unresolved, 
as markers of the pulsar ro tation rate. The p ulsar 
timing package TEMP02 ( Hobbs et al. 20061) was 
used to obtain a pulsar folding ephemeris for the 
observations. 


3.2. MWA data analysis 

3.2.1. Incoherent beamforming and de-dispersion 

Post-observation, the voltage samples for the 
two polarizations across the 128 tiles were squared 
to form powers and then summed in each 10 kHz 
channel for each 100 ps time step to form an 
incoherent beam. The sum was inverse-variance 
weighted to maximise the signal to noise. This in¬ 
coherent summing of the MWA tiles yields a large 
field-of-view, 24°(each tile beam oc A/d, where A is 
the observing wavelength, and d is the size of the 
tile), but a sensitivity that is \/]V times smaller 


than that achievable by coherent combination of 
the full array, where N is the number of tiles. 


The incoherent beam was then convert ed int o 

2004 . 


the PSRFITS data format (Hotan et al 


The files were then processed using the pulsar ex¬ 
ploration and search toolkit (PRESTO; Ransom 
2001). We averaged in time to 400 ps, de-dispersed 
the data, and generated a time series at the nom¬ 
inal DM of the Crab pulsar (56.70 pc cm -3 ). 



time since the beginning of the observation (s) 


Fig. 1. Dedispersed time series for the strongest 
giant pulse in the MWA data. The yellow (lighter) 
portions were used to estimate the off-pulse base¬ 
line shown by the red line. The gray (darker) area 
is the region that was integrated over to measure 
the total pulse energy. The plot spans 180 ms and 
has a time resolution of 400ps. 


Figure [I] shows the brightest giant pulse from 
our observations, where the scattering tail extends 
to ^50 ms, which is longer than the pulsar rota¬ 
tion period of 33 ms. Further analysis of scatter¬ 
ing and comparison with the previously published 
observations are discussed in § 13.41 


3.2.2. Measuring the properties of giant pulses 

The Crab Nebula is an extended radio source, 
measuring about 5.5' in diameter, with a flu x 
density of ~955i/ -0 27 Jy dBietenholz et al. 1997 ). 
where v is the frequency in GHz. In our case, 
the Nebula occupies only a tiny fraction (~0.2%) 
of the telescope beam. The nominal sky tem¬ 
perature at our observing frequency and for our 
pointing was evaluated from the global sky model 
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( de Oliveira-Costa et al.l I2008T ) to be 243 K; this 


value explicitly includes the contribution from the 
Crab Nebula. At 200 MHz, our receiver tempera¬ 
ture is 25 K (lTingav ct al . 20131) . The incoherent 
gain ( Gi nco ) of the instrument is given by the ex¬ 
pression 


Girt.ro — 


£( 16 Vn) 


2k 


2 0 , 


( 4 ) 


where A is the observing wavelength, N is the 
number of antenna tiles, k is the Boltzmann’s con¬ 
stant, and 6 is the elevation angle. For a nomi¬ 
nal Gi nC o of 0.012 K Jy -1 at the Crab’s elevation 
(41.42°), the system equivalent flux density (S sys 
= T sys /Gj nco ) of our telescope is ~22000 Jy. 

The phase of each giant pulse was determined 
from the arrival times of the giant pulses and the 
pulsar ephemeris derived from the Parkes data (see 
§ED- We find giant pulses only at the MP and 


earlier results (Lunderen et al. 

1991 Cordes et al. 

2004; Soerlasnov et al. 2004; 

Podov et al. 

200! 

Bhat et al. 20081: Karunnusamv et al. 2012 

). 


The total flux density was measured by inte¬ 
grating over a 40 ms window for every pulse pe¬ 
riod, beginning from the pulse arrival time. The 
integration window was determined by estimat¬ 
ing the timescale over which the amplitude of the 
brightest pulse fell below the 1 a noise level of 
the baseline. This corresponded to 40 ms, ap¬ 
proximately 6 e-folds of the estimated scattering 
timescale (see §& The baseline was determined 
by a linear least square fit to the off-pulse power 
across a series of five windows, each 12 ms wide, 
distributed before and after every pulse period in 
the timeseries (see Figure 1). 


The resulting distribution of integrated flux 
density is consistent with the noise distribution 
of the integrated timeseries, because the number 
of giant pulses are few compared to the number 
of the observed periods (~80,000 periods). The 
root-mean-square noise variation (rms) was deter¬ 
mined by estimating the standard deviation of the 
distribution, and the S/N was subsequently deter¬ 
mined. The integrated flux density of each giant 
pulse was then calibrated using equation 0 


In order to determine the number of false posi¬ 
tives our pipeline generated, we simulated a noise 
dominated timeseries using the mean, and the 
variance of our data. We then subjected the sim¬ 


ulated timeseries to the same analysis described. 
The number of false positives generated by this 
experiment above 3.5<r was two, a false positive 
rate we deemed acceptable where a is the rms of 
the distribution. At this threshold, our minimum 
detectable flux density is ~70 Jy. 


In the actual analysis of the observed data-set, 
above a threshold of 3.5cr, 45 giant pulses were re¬ 
covered from the MP phases and 10 giant pulses 
were recovered from the IP phase. The majority 
of giant pulses, 84%, occur at the phase of the 
MP and the remaining 16% occur at the IP phase. 
This is in agreement with previous resu lts pub¬ 
lished in the literature (Popov fc Stannersl 2007t 
Bhat ct al. 2008 : Mickaligcr et alJ 2012h . 


3.3. The giant pulse fluence distribution 


We now make the assumption, common in the 
literature, that the measured giant pulse energy is 
drawn from a power-law distribution, and attempt 
to determine the power-law index. We are in fact 
measuring the time integrated flux density, or flu¬ 
ence, which is a quantity directly proportional to 
pulse energy received by a given collecting area 
and will use this quantity in our analysis. 


A power-law index is often determined geomet¬ 
rically by a least squares fit of a strai ght line 
to logarithmically binned data ([Bhat et ah 20081 ). 
This is not the most reliable method of estimat¬ 
ing this parameter and is generally considered 
to generate biased values of the power-law index 
( Goldstein. Morris fc Yenl 12004 ). In order to re¬ 
move subjectivity and to limit bias we have chosen 
to use Hill’s estimator as a maximum likelihood 
estimator (MLE) of the index of the underlying 
fluence distribution. With n measurements of flu¬ 
ence, Xi , we determine the index, /?, of the under¬ 
lying power-law distribution via: 


b 


1 + n 



( 5 ) 


Hill’s estimator is sensitive to the determina¬ 
tion of the power-law cut-off (a; m in), but geometric 
methods are also very sensitive to this and often 
the power-law cut-off is chosen somewhat subjec¬ 
tively. To limit this subjectivity we have chosen 
to follow a method of minimising the Kolmogorov- 
Smirnov distance between the fitted distribution 
and the measured distribution, as a function of 
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the power-law cut-off, in order to determine the 
optimum value of a; m in- The standard error in the 
power-law index estimate, a = (/? — was 

also determined analytically. Using this method 
we determined the power-law index of fluence dis¬ 
tribution detected at the MWA to be —3.35 ± 
0.35. 


The cumulative distribution (CDF) of the 
main-pulse phase giant pulses detected with the 
MWA is shown in Figure [5J The slope of the CDF 
of giant pulse fluences can also be used to deter¬ 
mine the index of the underlying power-law distri¬ 
bution, although error prone, this method is often 
used in the pu lsar l iterat ure ( Popov fc Stappersi 


2007 : iKaruppusamv et ah 2Qlol h The gradient is 
equal to 1 + /3, but as the index is formed from a 
least squares fit of a straight line to logarithmic 
bins that are not independent, the error in the 
index is difficult to determine. In this case we find 
a gradient of —2.22 which would predict a /3 of 
—3.22, consistent with our MLE analysis. 

We applied the MLE method to our observa¬ 
tions of the Parkes pulses (at 1382 MHz) and ob¬ 
tained an index of —2.85 ± O.Off^ which is con¬ 


sistent with that obtained bv IKaruppusamv ct al 


(1201 Oil : it is slightly steeper than the —2.33 ± 
0.14 and 2 20 ± 0.18 a t 1300 and 1470 MHz re¬ 
ported by iBhat et HI ( 20081) . but is within the 
range of —2.1 ±0.3 t o —3.1 ± 0.2 as given by 
Mickaliger et al. ( 2012h . 


3.4. Pulse shape and scattering 


Scattered pulse shapes are typically modelled 
as convolutions of intrinsic pulse shapes with the 
broadening functions characterising multi-path 
scattering through the ISM and the instrumental 
response. By a pplying a d econvolution method 
as described in IBhat et al.l (200311 to the pulse 
shown in Figure [Q we estimate a pulse broaden¬ 
ing time, Td ~ 6.1 ± 1.5 ms, which can be com- 
pared to ~ 0.67 4= 0.10 ms reported previously by 


Bhat et al.l (j2007i ) from observations made with 
an early prototype built for the MWA. Our es¬ 
timated Td is thus nearly 10 x larger than that 
measured earlier with the 3-tile prototype MWA 
and 5 x larger tha n the extrapolated val ue fr om 
the observations of IKaruppusamv et al. ( 2012 1 at 
a nearby frequency of 174 MHz. However, it is a 
factor of two less than that expected based on re¬ 


cent observations bv lEllingson et al~l ( 2013 1 with 
the LWA and their re-derived freque ncy scaling 
i ndex Td o c t/~ 3 7 . As high l ighted bv IBhat et al.l 
( 2007 1 and Ellingson et al. I ( 2013ll . the Crab pul¬ 
sar is well known for highly variable scattering, on 
timescales of the order of months to years, which 
can be attributed to ionised clouds or filaments 
within the nebula crossing the line of sight. 


3.5. Coincidences in the MWA and Parkes 
giant pulse arrival times 


A search for coincident giant pulse arrival times 
was performed. Of the 55 giant pulses that were 
detected with the MWA, 45 were detected dur¬ 
ing the period of simultaneous observations with 
the Parkes radio telescope, and 1681 giant pulses 
were detected at Parkes. Due to the dispersive 
nature of the interstellar medium, the high fre¬ 
quency component of broadband giant pulses will 
arrive at the Parkes observing frequency earlier 
than at the MWA fr equency. The Jo drell Bank 
monthly ephemeri^ dLvne et al.lll993h provides a 
DM for the Crab pulsar of 56.774±0.005 pc cm -3 . 
The relative delay in pulse arrival time for this 
DM is therefore 5.746±0.008 s between 1382.00 
and 200.32 MHz (upper band edge for the MWA). 

In order to search for coincident giant pulses ob¬ 
served at the MWA and the Parkes radio telescope, 
we performed a correlation analysis between the 
arrival times of the MWA and the Parkes giant 
pulses. This analysis revealed that the backend 
used for data recording at Parkes had a clock er¬ 
ror of 36.000 s in the timestamps at the time of 
this observation. As a result, the timestamps for 
the Parkes pulses are ahead of UTC by 36.000 s. 
We found 23 coincident giant pulses at —30.248 s 
(see Figure O. Correcting for the 36.000 s error 
yields 5.752 s delay (MWA trailing Parkes), con¬ 
sistent with the dispersive time delay between the 
arrival times of the MWA and the Parkes giant 
pulses. 


The 23 coincident giant pulses detected at the 
MWA and Parkes implies 51% correlation, since 
45 giant pulses were detected with the MWA dur¬ 
ing the common observing period. Our observa¬ 
tions are ther efore consistent with previous obser¬ 
vation s (e. g. Cornelia et al.1 ( 1969h : Sallmen et al 
(Il999h l. see ? 14.21 for comparison. 


2 http: //www.jb. man. ac.uk/pulsar/crab.html 
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3.6. Spectral indices of coincident giant 
pulses 


The spectral indices of the giant pulses detected 
simultaneously at the MWA and the Parkes radio 
telescope were estimated using oc u a , where 
is the fluence of the giant pulse at frequency 
v and a is the spectral index. The values de¬ 
rived here are presented in Figure QJ and display 
a range between —3.6 and —4 .9, which are steeper 
than the limits obtained by iKaruppusamv et al.l 
(l2010h . -1.44 ± 3.3 and -0.6 ± 3.5 for the MP 
and IP, respectively. Our limits are consistent with 


Sallmen et al. ( 1999 1. based on their simultaneous 


detections of 29 giant pulses at 610 and 1400 MHz. 
Sallmen et al. ( 1999h constrain the spectral index 


spread to be within the range —2.2 > a > —4.9. 
The difference in the lower limit is discussed in 


4. DISCUSSION 


4.1. The fluence distribution 

4-1.1. Comparison with the literature 

We detected a total of 55 giant pulses at the 
MWA, and estimate a power-law index of >3 = 
—3.35 ± 0.35 for their pulse fluence distribution, 
using the maximum likelihood estimator method. 
For the Parkes data, in which we detected a total 
of 2075 pulses, application of this method yields 
(3 = —2.85 ± 0.05. In the published literature, 
however, the reported power-law index for the gi¬ 
ant pulse fluence distribution covers a wide range. 
Different distributions have been considered, in¬ 
cluding pulse peak amplitude, average flux density, 
and pulse fluence, although a degree of propor¬ 
tionality exists between these different measures 
direct comparison is difficult. Our power-law in¬ 
dex of —3.35 ± 0.35 for the pulse fluence distri¬ 
bution at 193 MHz is steeper than the range from 
— 1.51 ± 0.05 to —2.39 ±0.12 that was reported 
by Karuppusamv et al. ( 2012 ) at 110-180 MHz. 
Cordes et al. ( 2004! ) estimate —2.3 for the index 
of giant pulse peak amplitude distribution at the 
main-pulse phase; their observation was made at 
430 MHz and they note that the distribution only 
obeys a power-law for a limited range of pulse flux 
density. However our estimated (3 is in good agree- 
ment with a value of — 3.46 ± 0.04 determined by 
Lundarcn et al.l (119951 1 for the average flux den¬ 


sity distribution from their extensive observations 
at 800 MHz. 

We propose that the majority of the published 
variation in power-law indices arises from the un¬ 
derlying distribution not following a simple power 
law over a large range of fluences and also from 
inconsistencies in the subjective determination of 
the power-law cut-off. The MLE based method 
that we used here has the virtue of removing some 
subjectivity from this analysis. 


4-1.2. Effects of scintillation 

Observed radio pulses are often influen c ed by 


interstellar scintill ation (ILvne fe Rickett 
Cordes et al.l 120041 ). Under unfavorable condi¬ 


19681: 


tions, this effect can strongly modulate the ob¬ 
served fluence of giant pulses. In our case, this 
effect may arise from refractive modulation since 
short-term diffractive modulations will be highly 
quenched over our observing duration and band- 
width, at both M WA and Parkes frequencies. 
Cordes et al. ( 20041 ) report A v& < 0.8 MHz and 
Atd = 25 ± 5 secs at 1.48 GHz; where A Vd and 
A td are the diffractive scintillation bandwidth and 
timescale respectively. At the MWA frequency 
these will be orders of magnitude smaller and are 
therefore no longer relevant . As for refractive 
modulation, Rickett fc Lvne i 199 0; report a time 
scale of ~ 83 ± 50 days at 196 MHz and ^6—12 
days at their highest frequency of 610 MHz. Even 
if we consider a factor of two variability as a pos¬ 
sible worst case scenario at our frequencies, it 
may not still significantly alter our spectral index 
analysis and conclusions. 


4.2. Correlated detections and the intrin¬ 
sic giant pulse spectral index distri¬ 
bution 


Simultaneous observations of the Crab pulsar at 
multiple frequencies have been rep orted i n the l it- 
erature. This includes the work of Cornelia et al 


( 19691 ) from early pulsar observations, at 74 and 
111 MHz, to obs ervations that reach fre quencies as 
high as 8.9 GHz ( Mickaliger et al.ll2012l) . The level 
of correlation presented d oes vary. Observations 
bv lEllingson et al~l ( 20131 ) using the LWA with 4 
x 16 MHz observing bands spread from 28 MHz to 
76 MHz contain 33 giant pulse detections, only one 
of which was observed in all 4 bands, but almost all 
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were observed in at least 2. Based on observations 
with the Green B ank 25-metre telescop e and the 
Very Large Array. [Sallmen ct~aL (1993) reported 
the detection of 70% of giant pulses seen simulta¬ 


neous ly at both 600 and 1400 MHz. IPopov et al 


( 20061) report 27% simultaneous detections at their 
observing frequencies of 2 3 and 6 00 MHz bu t only 
16% at 111 and 600 MHz. Bhat et al. ( 20081 ) found 
70% correlation between the the giant pulses ob¬ 
served at 1300 MHz and 1470 MHz , simi lar to 
the result obtai ned by Sallmen et al. JH999T) . The 
observations of Mickaliger et, al.l ~( 2012h at larger 
frequency separations show a significantly smaller 
fraction, 3% to 5% between 1.2 GHz and 8.9 GHz 
for the MP and IP giant pulses, respectively. 

Our observations display a correlation of 51% 
between the pulses detected at Parkes and at the 
MWA, but the range of spectral indices displayed 
by the detected pulses (—3.6 to — 4.9) is consid¬ 
erabl y narrower than reported by ( Sallmen et al.l 
Il999l) . It is likely that the paucity of pulses with 
spectral indices shallower than —3.6 is due to the 
sensitivity limit of the MWA telescope and conse¬ 
quently that the intrinsic degree of correlation is 
likely to be much higher than 50%. We hypothe¬ 
sise that the apparent reduction of correlation is 
due to the wide range of spectral index displayed 
by the giant pulses, coupled with the sensitivity 
limits of the two telescopes. 

In order to investigate this further, a Monte 
Carlo simulation was conducted using the Parkes 
data as the known distribution of giant pulse flu- 
ence. Attempts were made to predict the number 
of pulses detected coincidentally at the MWA and 
their spectral index distribution using a two sam¬ 
ple Komolgorov-Smirnov (KS) test. 

We first determined the most likely intrinsic dis¬ 
tribution of gia nt pulse spect r al ind ices. As re¬ 
ported by both ISallmcn ct al.l ( 1999 ) and Karup- 
pusamy et al. (2010), the scatter in individual 
giant pulse spectral indices is large. The hypothe¬ 
sis that the intrinsic giant pulse spectral indices 
could have been drawn from a uniform distri¬ 
bution between —2.2 and —4.9 was immediately 
discounted; it reproduced the range of detected 
spectral indices (—3.6 to —4.9) however it over¬ 
estimated the number of coincident detections by 
a factor of 20. The broad distributions favoured by 
Karuppusamy et al. (2010) were also discounted 
by this analysis, as although they could reproduce 


the number of coincident detections, the distribu¬ 
tion of predicted spectral indices were not com¬ 
patible with the observations. A normal distri¬ 
bution of spectral indices was a much more suc¬ 
cessful hypothesis. It simultaneously satisfied the 
constraints of the detected spectral indices and 
the number of detections. We chose to investi¬ 
gate this further and determine the parameters of 
such a distribution that would best fit our data. 
We allowed the mean of the distribution to vary 
between —2.0 and —3.6, but held the standard 
deviation to 0.6, consis tent with the observations 
( 1999t) . and investigated the pre¬ 


of Sallmen et al 


dieted number of coincident detections and their 
spectral indices. 


For each of the Parkes giant-pulse detections we 
randomly drew a spectral index from the distribu¬ 
tion under test and predicted the fluence of the 
pulse at the MWA. If the predicted fluence was 
above the MWA detection threshold we registered 
a detection and the spectral index that was asso¬ 
ciated with it. This was repeated for every pulse. 
Each trial distribution was tested 1000 times and 
the results combined. 


As we have recorded observations at both tele¬ 
scopes we can reverse this analysis and use the 
MWA detected pulse and spectral index distribu¬ 
tion to predict the coincident pulses at Parkes, 
and thereby determine the percentage of corre¬ 
lated pulses. We are only observing the steep 
spectrum end of the giant pulse distribution at 
the MWA, as evidenced by inspecting the mea¬ 
sured spectral index distribution. We can infer 
that some pulses detected at the MWA will have 
a combination of fluence and steep spectral index 
that would bring them below the detection thresh¬ 
old at Parkes. This feature would be the reason for 
the less than 100% correlation in pulse detections. 

In Figured! the lower panel shows the predicted 
number of GPs detected at the MWA for each trial 
spectral index distribution as triangles and the 
percentage of coincident GPs between the MWA 
and Parkes. There are also two lines representing 
the observed number of GPs and the observed co¬ 
incidence fraction. To satisfy the observations, a 
trial spectral index distribution should simultane¬ 
ously satisfy both these constraints. The error bar 
is the standard deviation about this average pre¬ 
diction. The upper panel shows the results, both 
distance and p-value, for the two-sample KS test 











































between the predicted spectral index distribution 
of the coincident pulses and the measured distri¬ 
bution; the error bar is the standard error. The 
p-value can be interpreted as the probability that 
the distance is as large as observed, if the null hy¬ 
pothesis is true (being that both populations are 
drawn from the same distribution). The null hy¬ 
pothesis is typically not rejected if the distance is 
below a critical value for the two distributions un¬ 
der test. These results show that, for those distri¬ 
butions with a significant p-value, we cannot reject 
the null hypothesis and that the observed and sim¬ 
ulated spectral index distributions are consistent 
with same parent distribution. 

Without any fine tuning, the simulations pre¬ 
dict the following: 

• The detected spectral index distribution at 
the MWA is best replicated by an input spec¬ 
tral index distribution of -2.85 with a stan¬ 
dard deviation of 0.6 

• The observed 50% correlation is best repli¬ 
cated by a spectral index distribution closer 
to -3.15 

• Both these input distributions over-estimate 
the number of coincident pulses by a factor 
of 2 to 3 

4-2.1. Fine tuning the analysis 

The simplest, and most reasonable fine tuning 
is to reflect the lack of flux density calibration and 
assume some less-than-theoretical detection sensi¬ 
tivity. We can align the number of MWA detec¬ 
tions with the best fit spectral index distribution, 
using the MWA detection threshold as a fit pa¬ 
rameter. Raising the MWA detection threshold by 
10% generates a best fit mean spectral index distri¬ 
bution of -2.8, and correctly predicts the number 
of MWA detections. However it cannot simultane¬ 
ously replicate the correlation percentage, which is 
now reproduced by invoking a source spectral in¬ 
dex distribution with a mean of -3.2. 

Altering the width of the model spectral in¬ 
dex distribution does not provide the required fine 
tuning. The width of 1.2 is wide enough that it 
produces features similar to that generated by a 
uniform distribution of spectral indices. In that 
we can reproduce the number of detections at the 


MWA but the distribution of observed spectral in¬ 
dices differs too greatly. In the case of the narrow 
distribution, it is difficult to reproduce the num¬ 
ber of detections unless we adjust the relative sen¬ 
sitivity of the telescopes beyond what would be 
considered reasonable. 

A distribution that satisfies the constraints in 
predicting the Parkes pulse survival to the MWA 
observing frequency is therefore still inconsistent 
with the distribution that satisfies the degree 
of correlation. A more sophisticated fine-tuning 
would be to invoke some evolution in the distri¬ 
bution as a function of frequency. We have found 
that this evolution does not have to be extreme in 
order to remove the inconsistency. A 5% flatten¬ 
ing in the spectral index, on average, between the 
two frequencies is all that is required to match the 
observed pulse coincidence rate, while maintaining 
a spectral index distribution consistent with that 
observed. We obtained this estimate using the 
previous Monte Carlo analysis, but including an 
average reduction factor in the randomly drawn 
individual pulse spectral index when predicting 
the observed fluence at the MWA. We incorpo¬ 
rated the same factor in the inverse experiment. 

Incorporating this fine tuning, our observations 
are consistent with the hypothesis that the spec¬ 
tral index of giant pulses is drawn from normal 
distribution of standard deviation 0.6, but with 
a mean that displays an evolution with frequency 
from -3.00 at 1382 MHz, to -2.85 at 192 MHz. 
Note this estimation only includes the integrated 
effect of the evolution in spectral index, and there¬ 
fore does not incorporate curvature. It would be 
consistent with these results if the low frequency 
spectral index were considerably flatter, or the 
high frequency steeper, than that estimated here, 
but that would require a rapid evolution, such as 
a turn-over, to be present. Figure [B] shows the re¬ 
sults of the analysis including the fine tuning and 
we summarise the results here: 

• The MWA detection threshold is 10% less 
sensitive than our theoretical prediction 

• There is a 5% flattening in the mean of the 
spectral index distribution from -3.00 to - 
2.85 between 1382 MHz and 192 MHz, the 
spectral index distribution in this analysis 
has width of 0.6, but this has not been con¬ 
strained by the fits 
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• This distribution correctly predicts the num¬ 
ber of coincident detections and the 50% cor¬ 
relation between the two observations 


The normal emission from pulsar s typic ally 
obey s a simple power law of index —1.8 (jMaron c t al 
200 0). However, the MP and IP components of 
the Crab pulsar have much steeper spectral in¬ 
dices of -2.8 and -3 .7 when measured bet ween 
410 and 1660 MHz ( Manchester et al. 197 lh . A 
flattening in the spectra of normal pulsar emis¬ 
sion is commonly observed and many young pul¬ 
sars have been found to display complex spectra 
and low frequency (~100 MHz) spectral tur n-overs 


(jMalofeev et al.l 11994 iKiiak et al. 120071) . The 
phase averaged emission from the Crab pulsar is 
complicated due to the multi-component struc¬ 
ture. We are unaware of any observed flattening 
of MP component at lower frequencies but sus¬ 
pect that this is confused by the rapid evolution 
of precursor. The cause of the spectral turnover 
in these objects is not known and experiments are 
complicated by the difficulty of observing complex 
pulse profiles at low frequencies. 


4.3. Propagation effects within the nebula 

Our measurements of the scattering timescale 
( 43.41) indicate that the environment of the neb¬ 
ula continues to affect the observations of the p ul¬ 


sar, as a lready observed by [R ankin fe Counselman 
(1973), Isaacman.&Rankm (197 71), Bhat et al 


( 2008 ). Kuz’min et al.l ( 20081) and Ellingson c t al. 


(120131 ). Free-free absorption within the filamen¬ 
tary structures, evident in observations of the 
Crab nebul a, has already been repo rted byE n d 2 
( 1970h and Bietenholz et al.l ( 1997 ). We now in¬ 
vestigate whether this effect could be responsible 
for the modest flattening in spectral index that 
has been indicated by our analysis. The free-free 
absorption coefficient (a®) in the Rayleigh-Jeans 
regime can be numerically approximated (in cgs 
units) to be: 


a® = 0.018 T- 3/2 Z 2 n e mn- 2 
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( 6 ) 


where Z is the atomic number of the absorber, 
n e and n- u are the number density per cm - 3 of 
the electrons and ions respectively, and T is the 
temperature in K ( Rvbicki fc Lightmanlll979f ). gs 
is the Gaunt factor, which in this regime (the high 


te mperature or quantum limit ) is approximated 
bv I Cavetl ( 1970t ): 


gs = 


A /4fcT\ 

7 r n \hve"t) ’ 


(7) 


where fcis Boltzmann’s constant, h is Planck’s con¬ 
stant, and 7 is Euler’s constant. At these temper¬ 
atures and frequencies the Gaunt factor is approx¬ 
imately 10. The optical depth 7 # is given by: 



( 8 ) 


where z is the path length, and the absorption 

_ ff 

factor is subsequently calculated as e~ T . 

The nebula has an approximate radius of 1.7 pc 
and if all the observed DM (56.70 cm“ 3 pc) were 
due to the nebula, then n e would be 56.7/1.7 or 
33.4 cm -3 . Assuming the nebula is predominantly 
Hydrogen, Z = 1, n e = rq = 33.4 cm -3 at a tem¬ 
perature 10 4 K, then at our observing frequency 
(192 x 10 6 Hz) the free-free absorption coefficient 
is ~ 5.4 x 10 ~ 23 cm -1 . This is insufficient to ac¬ 
count for any absorption even at these very high 
Gaunt factors. 


The ionised material is not evenly distributed, 
and this inhomogeneity will strongly influence the 
level of free-free absorption. The filamentary 
structures of the nebula that have been observed 
to displa y free-free absorption in spectral index 
maps by Bietenholz et alJ ( 19971) . have been esti¬ 
mated to have an electron density in excess of 200 
cm -3 . These are arcsecond size structures, imply¬ 
ing a size scale of 0.02 pc. Each of these filaments 
would contribute 1.2 x 10 -4 to the optical depth 
and attenuate the intensity by a factor of 0.9998, 
which also cannot supply the necessary absorption 
even if s everal hundred filaments in tercept the line 
of sight. iDavidson fc Tucker (19701) use line inten¬ 
sities to estimate the filamentary n e to be ~ 1000 
cm -3 , however each filament now contributes ~ 10 
cm~ 3 pc to the dispersion measure so only a small 
number of crossings are permitted, therefore these 
structures cannot provide the level of absorption 
required. 

There is also evidence of very small (AU) scale 
filamentary struct ures with electron de nsities in 
excess of 10 4 cm -3 ( Smith fe Terrvll201ll ). If many 
tens of these small filaments are intercepted, they 
are capable of generating the few percent level ab¬ 
sorption required (few %). They would also not 
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contribute more than 10 cm~ 3 pc to the pulsar dis¬ 
persion measure. However it is unlikely that these 
structures ex ist in sufficient numbe rs to contribute 
at this level ( Smith fe Terrvlfeoil ). 


In summary it is unlikely that free-free absorp¬ 
tion within the filamentary nebula is capable of 
attenuating the pulse intensity sufficiently to ex¬ 
plain our predictions, and therefore any variation 
in the spectra of the giant pulses is unlikely to be 
propagative in nature. 


5. CONCLUSIONS 


In this article, we report observations per¬ 
formed simultaneously with the MWA, at a fre¬ 
quency of 193 MHz (15 MHz bandwidth) and the 
Parkes radio telescope, at a frequency of 1382 
MHz (340 MHz bandwidth). We detected a to¬ 
tal of 55 giant pulses at the MWA, and estimate 
a power-law index of /? = —3.35 ± 0.35 for their 
pulse fluence distribution. 23 of the pulses were 
detected at both the MWA and at Parkes ra¬ 
dio telescope. These results are consistent with 
the spectral index of the giant pulses being drawn 
from normal distribution, the parameters of which 
are frequency dependent. The mean of the spec¬ 
tral index distribution varies from —3.00 to —2.85 
between 1382 MHz and 192 MHz. It is unlikely 
that this flattening can be caused by any propaga¬ 
tive effects within the nebula. 


In this work we have proposed that the less 
than total correlation observed is simply a func¬ 
tion of the spread in giant pulse spectral indices 
and the relative sensitivities of the two instru¬ 
ments, but have not proposed a mechanism by 
which this large range of spectral indices may 
be generated. Giant pulses are thought to com¬ 
prise complex temporal substruct ures on nano¬ 


secon d to micro-second ti mescales (jHankins et al 
2003h . Eilek et al. ( 2002 1 have noted that these 


substructures are strongly correlated at two ob¬ 
serving frequencies with a small fractional band¬ 
width, and less correlated as the fractional band¬ 
width increases, suggesting a complex behaviour 
as a function of frequency, cannot be ruled out by 
the behaviour we see in this experiment. 


The MWA’s large frequency coverage and flex¬ 
ible system design offer unique opportunities to 
further investigate the wide-band properties of gi¬ 
ant pulses. With the full-bandwidth VCS record¬ 


ing now feasible (ITremblav et al.N2015l) . it is pos¬ 
sible to perform sensitive observations simultane¬ 
ously at multiple frequencies within the 80 - 300 
MHz range, by suitably spreading out the 30.72 
MHz observing band across a large range in fre¬ 
quency. This will allow us to conduct observa¬ 
tions that span large fractional band widths, with 
the prospects of determining the emission band¬ 
width and spectral nature of giant pulses at low 
frequencies. 


Combining all the MWA tiles coherently, we 
will realise a factor of ~11 improvement in sen¬ 
sitivity over this experiment. With an increase in 
gain by this factor, as well as increase in band¬ 
width by a factor of 2, the instrument is ex¬ 
pected to detect about 10,000 pulses in a one 
hour observation. This also implies that about 
600 pulses are expected to be detectable in a 1.28 
MHz subband of the MWA. This increase in sen¬ 
sitivity means that a higher percentage of coinci¬ 
dent pulses will be detectable when observing in 
the coherent mode with the MWA full bandwidth. 
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Fig. 2. Top: The cumulative distribution of the 
MP giant pulses detected with the MWA. Also 
shown here is MLE power-law fit to the distribu¬ 
tion, with a power-law index of —2.35 ± 0.35. The 
shaded region is a ±0.35 error in the power-law in¬ 
dex. Bottom: Distribution of the MP giant pulses 
detected with the Parkes radio telescope. The 
plotted errors of fit might be underestimated. The 
solid line is the MLE power-law fit to the data with 
a power-law index of —2.85 ± 0.05. The shaded is 
region is a ±0.05 error in the power-law index. 
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time delay (s) 

Fig. 3.- A plot showing the search for coinci¬ 
dent giant pulses in the MWA and the Parkes 
arrival times. We identify a 36.000 s clock er¬ 
ror in the timestamps of the recording machine 
used at Parkes at the time of this observation (see 
text). We used the arrival times of the Parkes gi¬ 
ant pulses to search for coincidence in the MWA 
pulse arrival times. As can be seen in the figure, 
we found strong peak corresponding to 23 pulses at 
—30.248 s. Correcting for the 36 s offset yields 5.75 
s which corresponds to the time delay between the 
arrival times of the MWA and the Parkes pulses. 



i.oo-| 


20 . 0 -| 


14 
























Monte Carlo Analysis Results 




Mean of Spectral Index Distribution 


Fig. 5.— The results of the Monte Carlo analy¬ 
sis to determine the best fit intrinsic mean spec¬ 
tral index of the giant pulses. The three con¬ 
straints, number of coincident detections (23), ob¬ 
served correlation factor (50%) and distribution 
of observed spectral indices cannot be satisfied by 
the same input distribution. 


Fig. 4.— The fluence in Jy-s of the coincident 
giant pulses detected with the MWA and Parkes. 
The error bars indicate the uncertainty in the sky 
temperature measurement for the MWA (which 
translate to ±0.08 Jy s in this case), which is 
negligible at the Parkes frequency. The vertical 
and the horizontal dash lines are the detection 
thresholds for the MWA and the Parkes obser¬ 
vations. The solid lines are spectral indices a = 
-3.6 an d a = -4.9. This figu re is similar to Fig¬ 
ure 7 of Sallmen et ah . •: 19!)l): i, in which they con¬ 
strained the spectral indices within the range of 
-2.2 > a > -4.9. 



Fig. 6.— Fine tuning the Monte Carlo analysis by 
incorporating an overestimate of the MWA sensi¬ 
tivity by 10% and a mild (5%) evolution in the 
spectral index as a function of frequency produces 
an internally consistent estimate for the mean 
of the spectral index distribution(—3.00 at 1382 
MHz, to greater than —2.85 at 192 MHz). 
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